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Abstract 


Using a fully self-consistent quantum statistical model, we demonstrate the 
possibility of Bose condensation of alpha particles with a concomitant phase 
transition in heavy ion collisions. Suggestions for the experimental observation of 
the signature of the onset of this phenomenon are made. 
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Theory and Results 

For some time, there has been an intrinsic interest in Bose condensation of 
particles with a concomitant phase transition in nuclear physics. Studies of 
multifragment final states over a broad range of incident energies may provide 
information about liquid-vapor and other types of nuclear phase transitions, and 
ultimately information about the equation of state of nuclear matter (refs. 1 through 
12), provided that multifragment disintegrations are related to bulk instabilities of 
nuclear matter at low densities. Most experimental observations (ref. 8) provide end- 
point information and do not necessarily provide information about transient 
phenomena. Here we demonstrate that Bose condensation of alpha particles is such 
a transient phenomena, which has thus far eluded experimental observation. The 
aim of this letter is two-fold: first, to establish the possibility of Bose condensation of 
alpha particles and the concomitant phase transition; second, to suggest specific 
experimental observations for its verification. The smallness of the finite volumes 
being considered could mask a clear signal of the phase transition. However, by 
exercising sufficient care, experimental signatures for the onset of alpha 
condensation may be observable. 

Bose condensation of pions was predicted earlier (ref. 6) to occur in the 
expansion phases of nuclear collisions. There is, however, a subtle difference 
between Bose condensation of alphas and that of pions. In the formation of alpha 
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particles, one requires four nucleons to be present in close physical and momentum 
proximity. Therefore, matter present must have some minimum density for 
condensation to occur. On the other hand, pion production is basically a one- 
nucleon phenomenon and therefore does not depend upon the presence of a 
minimum density of nuclear matter. Hence, pion condensation is observable in 
end-point measurements, such as total yields or cross sections, but alpha 
condensation is not. Its observation must be related to pre-ultimate stages of 
rarefaction (e.g., compression). 

In this letter, we show that for a given excitation energy (represented by a 
temperature T), one must reach a critical density (p c ) to observe the phase transition. 
Although the condensed alphas are present at densities below p c , their presence is 
masked by other (noncondensed) alphas produced in the collision. Our results 
indicate that the number of condensed alpha particles increases as the density 
increases toward p c . Therefore, our approach to investigating this phenomenon is 
to start at minimum (zero) density and compress the nuclear matter in the collision 
until the critical density (where the system becomes unstable) is reached. 

To focus on the physics of the phenomenon, a simple model is assumed. 

More detailed and reliable predictions would certainly require a more sophisticated 
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and realistic model; nevertheless, a simple model will suffice to demonstrate the 
phenomenon. 

We begin with an initial assembly of nucleons consisting of Nj n neutrons and 
Zi n protons which have been compressed and are evolving in thermodynamic 
equilibrium at an excitation energy (temperature) T to form clusters of light ions 
( 4 He - a, 3 He - h, tritons - 1, deuterons - d) together with neutrons (n) and protons 
(p). (Although the existence of a fully equilibrated system in heavy ion collisions is 
still debated, it is generally believed (refs. 4, 5, 9, 10, 11, 12) that equilibrium 
thermodynamics is a viable theory which gives reliable predictions). In the model, 

charge and baryon number are conserved using 

N p + N d + 2N h + N t + 2N a = Z m 

N n + N d + N h + 2N t + 2N a = N*, (1) 

where the Ni refer to the number of particles of species i. These particles have finite 
dimensions and initially move within a volume Vi n . For convenience, we treat the 
particles as point particles (no physical dimensions) which move in a reduced 
volume given by 

v pc = v »-X N i v i, ( 2 ) 

i 

where Vi is the eigen volume of the i— particle. The number (matter) density is 


P = (Nin + Zin)/V in - 


(3) 



The density of point particles is 
Ppt = (Njn + Z|n)/Vpt , 
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( 4 ) 

where p pt > p since V p t < Vi n . Our results will be presented in terms of p pt . 

For a distribution of fermions, quantum statistics (ref. 13) yields 

N i = 2g,V pt)t - 1 /2>.; 3 F FD (n/k B T), (5) 

where gj = 2 Si + 1 is the spin degeneracy, Pi is the chemical potential, h is the 

thermal wavelength for the i th particle of mass m, 

X i = 2 7ch(2jcm i k B T)- 1 /2, (6) 

where k B is Boltzmann's constant, and the species index is i (= p. n. h, or t) as 

appropriate. Values for the Fermi functions (spin 1) 

(7) 

Ffd (v) = I — ^ — dx, 
y ’ I 1 +e x v 

Jo 

are tabulated in the literature (ref. 14) and are increasing functions of v. For bosons 
(refs. 13, 15, 16), we have 

Ni = [exp (-pi/k B T) - 1] 1 

+ gi V pt Xj 3 F B e (-pi/k B T), 

where the first term on the right side of equation (8) gives the number of condensed 
bosons and the species index is i (= d, a). 
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In equation (8), the Bose-Einstein functions are given by 

(9) 

Fbe M = ! dx ‘ 

In equilibrium, the chemical potential for the i lh species is 

Pi = Zi pp + iij p n + e i, (10) 

where z, and nj are the number of protons and neutrons in the i th species and e * is 
its binding energy. This completes the description of the model used in this work. 

To illustrate the application of these methods to a model system, consider an 
assembly of Nj n = 40 neutrons and Zj n = 40 protons, such as might exist in a 
^Ca - ^Ca central collision. Figure 1 displays the distributions of light clusters for T 
= 1, 2.5, 5, and 7.5 MeV as a function of the density ratio p p t/Po where p 0 = 0.17 fm' 3 . 
Note that the distribution of clusters is completely dominated by alpha particles at 
these excitations. Other particles /clusters begin to appear at higher temperatures 
and lower densities. From the figure, we observe that the total number of alphas 
(condensed and other) remains constant as the system is compressed, i.e., as density 
is increased. However, the percentage of condensed alphas increases with increasing 
density (i.e., compression) so that at higher densities (for a given temperature) the 
dominant contribution comes from condensed alphas. It is these condensed alphas 
which ultimately control the highest stable density (critical density p c ) beyond which 
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instability sets in. This instability is evident in our results by a lack of stable 
solutions for densities greater than the critical values. The critical density ratios 
(Pc/po) for the onset of instability at T = 1, 2.5, 5, and 7.5 MeV are 0.5, 2.1, 6.1, and 11.3, 
respectively. For comparison purposes, we note that the corresponding freeze-out 
density ratios (taking into account the finite cluster volumes) are 0.5, 2.05, 4.61, and 
4.29, respectively. 

Mathematically, the instability is related to the singularity resulting from 
p — > 0 in the first term on the right-hand side of equation (8). Physically, the 
instability is related to the phase transition, which would be observable if the 
nuclear systems were large enough (much greater number of particles in the 
ensemble), and if the system could be observed for a long time. In ideal quantum 
gases, this phase transition occurs with Bose condensation. At this point, alpha 
particles could be created with little or no energy investment. In our work, it is 
important to recognize that the alpha condensation occurring during compression is 
a hot one. Its existence is not limited to near zero temperatures. A remarkable 
feature is that this condensation, if it occurs, should be experimentally observable. 

The question is, how do we observe it? The most common observable for any 
distribution of particles is to measure total yields /cross sections. A measurement of 
the total yield of alphas, however, would not confirm or deny the existence of 
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condensed alphas. From Figure 1, we note that the total yield of alphas (condensed 
and noncondensed) is constant and essentially independent of the density when 
below the critical density. Only the relative proportions of condensed to 
noncondensed alphas vary. Consequently, experimental observation/ detection 
must be based upon methods which distinguish between these two categories of 
alphas. One such method of detection might involve measuring the momentum 
distributions of the alphas. In the center-of-mass system, the condensed alpha 
particles will be produced with near-zero momenta, whereas the noncondensed 
alphas will have much larger momentum values. Therefore, a possible signature 
for alpha condensation might be found in experimental measurements of total 
transverse momentum (P ± ) or total in-plane transverse momentum ( P x ) of alphas 
which should decrease with increasing density. Ideally, for an infinitely-large 
system, the presence of a phase transition would be evident by a sudden 
(discontinuous) decline in measured momenta as the critical density is reached. For 
real, finite systems, any discontinuous reductions in measured momenta may be 
obscured by the small sizes of the systems considered; nevertheless, significant 
decreases in measured momentum components of alpha particles as the matter is 
compressed will provide indisputable evidence for Bose condensation of alpha 
particles with the concomitant phase transition. 



10 


Another experimental observable which may be linked to the Bose 
condensation of alpha particles is the possible occurrence of a critical maximum 
density p c for a given excitation (T), a density above which stable nuclear matter is 
non-existent. Recall from our model that the critical density arises mathematically 
from the singularity occurring in the term [equation (8)] representing Bose 
condensation. This is less definitive, however, than directly observing alphas with 
near-zero momentum as discussed earlier. 

Figure 2 displays the number of condensed alpha particles as a function of 
Ppt/Po for different temperatures. From the figure, it is apparent that the critical 
density where the instability occurs does depend on the temperature. Clearly, 
however, the ratio of condensed to noncondensed alphas at the critical density is 
independent of temperature. It is also interesting to note that for temperatures 
below the binding energy per nucleon of the alpha particle the number of condensed 
alphas reaches the same fixed value ( = 17 here) before the critical density is reached, 
irrespective of the system temperature. 

Finally, we emphasize that the findings and conclusions of the present work, 
although definitive and quantitative in nature, need further study and refinement 
using more sophisticated models than the one used here. Nevertheless, 
experimental verification of our proposal should be possible even though difficult 
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because of the finite sizes of nuclear ensembles. In view of the importance of the 
phenomena, verification of its existence merits concerted efforts by experimental 


investigators. 
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Figure 1. Distribution of particles up to critical density ratio for various temperatures. 

Curves for total alphas (condensed and noncondensed), condensed alphas and 
the next-most abundant other species are displayed. The initial ensemble 
contained 40 protons and 40 neutrons. 
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Figure 2. Number of condensed alphas as a function of density ratio at temperatures of 
1, 2.5, 5 and 7.5 MeV. Each curve terminates at the critical density ratio for 
that temperature. Note that the number of condensed alphas reaches the 
limiting value of 16.7 at the critical density, irrespective of temperature. 
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